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1.0 INTRODUCTION 


1 .1 Scope and Purpose of Report 

This progress report covers activities during the fourth 
quarter, March 25, 1976 through May 25, 1976, for LANDSAT 
Investigation #23790. This investigation is funded for 19 
months to develop techniques in Texas state agencies for using 
LANDSAT data to inventory and monitor coastal resources and uses. 
The General Land Office (GLO) is the Texas agency coordinating 
this investigation. Other participating agencies are the Bureau 
of Economic Geology (BEG), Texas Water Development Board (TWDB), 
and Texas Parks and Wildlife Department (TPWD). 

1 .2 Summary of VIork Performed 

During the fourth reporting period, activities centered 
around completing at least one scene in test sites 2 and 5 in 
humid and semi-arid zones of the Texas coast, respectively. In 
addition, plans were completed for the summer "test" of using 
image interpretation and ADP techniques to address General Land 
Office coastal management concerns, and for collecting data for 
the cost-savings analysis. 

Specific accomplishments include: 1) detailed analysis of 

image interpretation results from test sites 2 and 5, 2). refine- 
ment of the land use/land cover classification scheme developed 
for this investigation and correlation of this classification with 
the U.S.G.S. system, 3) introduction of a new program (H6R0UP) to 
the ADP classification schedule, and 4) preliminary development 
of a program to extract boundaries from ADP classification results. 
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2.0 PROBLEMS 

Anticipated examination of the effect of seasonal and other changes 
in 2 lANDSAT scenes for each test site (2, 3, and 5) has been delayed 
until late Summer or Fall because: 1) additional scenes for sites 2 

and 5 have not been completed due to anticipated additional technique 
development by June 1, 1976, when the analysis of site 4 1s scheduled 
to begin, and 2) there was a delay In receiving field support data for 
sites 2 and 5 from Texas Parks and Wildlife Department. 

3.0 ACCOMPLISHMENTS 

3.1 Data Acquisition 

The current status of LANDSAT data acquisition to support this 
investigation 1s reflected In Appendix A. The CCT's ordered for 
Test Sites 4 and 5 were received. Complete packages of imagery and 
digital data are now available as follows: Test Site 2 : 3 scenes 

(1 summer, 1 winter, 1 spring); Test Site 3 : 4 scenes (1 summer, 

2 winter, 1 spring); Test Site 4 : 4 scenes (1 suimier, 3 winter); 
Test Site 5 : 3 scenes (1 summer, 2 winter). As Indicated, these 

data sets provide some seasonal variations within each test site. 

In addition, for Test Sites 3, 4, and 5, the data sets within the 
winter season are separated by one to four years so that longer term 
trends In land use and land cover can be analyzed. 
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3.2 Examination of Test Site 2 (West Galveston Bay Area) 

3.2.1 Image Interpretation Results of Site 2 

Image interpretation of LANDSAT film transparencies reveals 
the complex active and abandoned natural drainage patterns, as well 
as dredged channels and the details of natural and jettied inlets 
within this coastal segment. San Luis Pass (Figure 1), an un- 
modified tidal inlet, is spanned by a highway bridge which is de- 
tectable with LANDSAT data. Extensive topographically low marshes 
just inside the pass were correctly identified and contain an 
abundance of Spartina alterni flora . Islands of less water content 
and higher reflectance vegetation surrounded by the low marsh 
were classified as high marsh, a decision which appears valid 
after examination of 1:40,000 color-infrared aerial photography 
(NASA Mission 325.) 

However, some initial errors in classification also occurred, 
and three examples are illustrated by Figure 1. A barely emergent 
marshy island (A, Figure 1) was mapped as highly turbid-shallow 
water. Where beach width decreases below about 80 m. (B, Figure 1) 
the beach, breaking waves, and the nearshore zone of turbid water 
appear as a zone of highly turbid or shallow water adjacent to 
the vegetated barrier flat. A strip of high-reflectance sand 
(C, Figure 1) is mapped as undifferentiated barren rather than 
dredge spoil because tht=» adjacent channel, about 50m. in w^dth, 
cannot be distinguished from the surrounding wet low marsh. 








The Freeport area (Figure 2) includes a major industrial -port 
complex nearly 3 km. in length and width, the urban areas of Free- 
port (A, Figure 2) and Surfside (B, Figure 2), and urban strip 
development (C, Figure 2) along State Route 332. This road leads 
northwest from Surfside but is not detectable where it crosses 
the low marshes landward of that city until the roadside develop- 
ment begins. At locality D (Figure 2), development becomes less 
continuous, but a string of individual industrial sites is a 
clue to the highway's location. On LANDSAT false-color composite 
transparencies, industrial sites are recognizable by: (1) the 

high-reflectance white to bluish-white tones caused by metal 
structures and the use of shell and sand fill, (2) the presence 
of holding ponds for liquids, and (3) distinguishable roads or 
dredged channels which lead to the site. Use of the first criteria 
alone can be misleading; at E (Figure 2) a rectangular area covered 
by barren dredge spoil approximately 160 by 400 m. fn size was 
misinterpreted as an industrial location. 

A detectable change in marsh type was mapped across a road 
(F, Figure 2) which, when checked in the field, was found to follow 
the crest of a hurricane-surge protection dike at an elevation of 
4.9 m. (16 ft.). Natural species zonation within saline marshes 
is related to gradual elevation changes, with less salt-tolerant 
species at greater elevations above mean sea level. In this 
instance, the dike (F, Figure 2) separates saline from brackish- 
to fresh-water marshes over an extremely short horizontal distance. 
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Figure 2. Freeport area, test site 2, with water 

units unclassified except for transportation 
(dredged canals or maintained river channel), 
delineated from LANDSAT scene 1289-16261, 

8 May 1973. Lettered locations are explained 
in text. 
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Spartina alterm' flora and Batis maritlma are abundant seaward of 


the dike, along with a relatively greater area of open water, while 
Spartina spartinae , Distichlis spicata , and some Scirpus spp. 
(bulrush) are found on the landward side. The false-color composite 
response of the latter area is dull red to blackish-red while the 
former is bluish black with faint reddish-black patches. 

3.2.2 Preliminary ADP Classification of Site 2 

Preliminary registered classification maps for the 8 May 
1973 scene (1289-16261) of Test Site 2 were completed during this 
reporting period. Because of the size and shape of the test area 
(7 quadrangle sheets arranged stepwise along the coast), each of the 
DAM, ELLTAB, CLASSIFY, and REGISTRATION steps required four runs 
to complete all 7 maps. This scene contains an exceptionally wide 
range of spectral levels, including particularly high reflectance 
areas, compared to those scenes previously studied. 

To improve the classification on this scene, the sample size 
for establishing the ISOCLS statistics was increased in some areas. 
After considerable experimentation, the parameters used were: 

1) maximum number of classes: 40, 2) maximum number of points 
per cluster;20, 3) maximum standard deviation:3.0, 4) minimum 
cluster separation: 2.0, 5) maximum iterations :10. After 
refinement of several clusters to reduce the standard deviations 
to acceptable limits,- 40 subclasses were retained. Due to limita- 
tions in the allowable size of the ELLTAB Table, the classification 
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was handled in tow groups of 25 and 15 subclasses respectively. 

The two groups will be consolidated manually to arrive at a single 
classification map for correlation with image interpretation results. 

After initial comparisons of the classification map with photo 
graphy (NASA Mission 300) and in conjunction with results of the 
HGROUP analysis (Section 3.5), the number of classes was reduced 
to 17. 

3.3 Examination of Test Site 5 (Southern Laguna Madre Area) 

3.3.1 Preliminary Description for Site 5 

In the vicinity of South Padre Island a pronounced deficit in 
precipitation leads to active aeolian transport and only sparse 
barrier island vegetation. Blowouts are common and storm channels 
bring wash-over deposits to Laguna Madre. Broad tidal flats have 
been formed by wind transport of sediment derived from the over- 
wash deposits (McGowan and Scott, 1975). Algal mats are present 
on parts of the wind-tide dominated flats; subaqueous vegetation 
(seagrass beds] occur within the lagoon. Marshes similar to those 
found in Sites 2 and 3 are scarce. Low areas along the lagoon 
margin support saline grasslands consisting primarily of Batis maritima, 
Monanthochloe littoral is and Borrichia frutescens (sea oxeye), while 
at greater elevations and in less saline soils (Johnston, 1955) 

Spartina spartinae is found. The fluvial -deltaic sands and muds 
of the Modern-Holocene Rio Grande delta system supports the culti 
vation of citrus, vegetables, grain sorghum and other crops. The 
Brownsville Ship Channel leads to shipping facilities and shrimp 
docks at the port of Brovins ville. 
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Test Site 5 shows greater complexity (Figure 3) than the 
other test sites--especially evident in the inactive Modern-Hol- 
ocene delta plains between the Rio Grande and the Brownsville Ship 
Channel and north of the Arroyo Colorado. The latter area is a 
complex of channels, tidal flats, subaqueous grass flats, algal mats, 
and undifferentiated barren substrate. The rectilinear boundaries 
in the region south of the Arroyo Colorado define cropland, seen 
on LANDSAT data as the gray of barren muddy fields, set within 
the dull brownish red of saline grasslands and the white to blue- 
white of undifferentiated barren clay-sand substrate. 

The city of Port Isabel (A, Figure 4) is detectable with 
LANDSAT data, as are the two causeways leading across Laguna Madre 
to South Padre Island. The lawns and trees of the city inter- 
spersed with paved roads and urban structures give a pebbly, dull 
red and white pattern which, especially when seen with transpor- 
tation features, helps identify urban complexes. Urban develop- 
ment on Padre Island, however, is not detectable because: (1) many 

streets are simply cleared areas of loose sand. (2) lawns and 
planted trees characteristic of developed areas are not present, 
and (3) structures immediately adjacent to the beach and within 
the foredune area are masked by the high reflectance of the barren 
sands. These residential and commercial areas therefore resemble 
the surrounding natural environment of barrier flat vegetation, 
tidal flat, and bare sand. A housing development consisting of 
dredged channels and spoil islands is mapped as tidal flat (B, 

Figure 4) since both environments appear simply as wet sandy sub- 

i 
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Figure 4. Brazos Santiago Pass area, test site 5, with 
classified land and water units (scene 2034- 
16205), lettered features are explained in text 



strata when individual channels cannot be detected. 

Also within the Brazos Santiago Pass area are subaqueous 
grass flats (C, Figure 4), barren sandflats (D, Figure 4) which 
are occassionally flooded (Groat et al., in progress), and saline 
grasslands. The long transportation feature, with associated dredge 
spoil, is the Brownsville Ship Channel. Note that the moderately 
turbid water issuing from the tidal pass turns northward since the 
nearshore circulation was under the influence of 12- to 16-knot 
winds from the south (U.S. Department of Commerce, 1975) at the 
time of satellite passage. 

As a result of work on this test site the land use/1 and cover 
classification scheme required further modification. The system 
shown in Table 1 is applicable to the entire Texas Coast and will 
be used during the remainder of the investigation. The numbers 
on the classification scheme refer to the revised U,S,G,S, system 
by Anderson, et al. (1976). 

3.3.3 Preliminary ADP Classification of Site 5 

Preliminary registered classification maps for the 25 
February 1975 scene (2034-16205) of Test Site 5 were completed. 

The area included five 7 1/2 minute quad sheets which were handled 
as three separate runs. Using the same statistics as for Site 2 
(Seciion 3.2.2), 25 classes were obtained. Based on the HGROUP 
analysis (Section 3.5), these were reduced to 17 classes. Maps 
showing both groupings were produced to be verified with aerial 
photography and maps of the Brownsvi 11 e-Harlingen Area (Groat, 
et al., in progress), and also in later correlation with the 
image interpretation results. 
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Table 1. 


LAND COVER AND LAND USE CLASSIFICATION 
FOR USE WITH UNENHANCED LANDSAT DATA 
IN THE TEXAS COASTAL ZONE 


1 Urban or Built-up Land 

2 Agricultural Land 

3 Grassland/Rangeland 

4 Forest Land 

5 Water 

(Level II omitted) 

6 Wetland 

7 Barren Land 


U-Mixed urban (16) 

Ui- Industrial (13) 

Ut-Transportation (14) 
Ue-Extractive-hydrocarbons (131 ) 

A-Cropland (21) 

G-Range-pasture (31) 

Gd-Vegetated dunes (311) 

Gb-Vegetated barrier flat (312) 
Gbr-Brushland (32) 

WO-Woodland or chaparral (43) 

WA-Non- turbid (501) 

WAst-Slightly turbid (502) 
WAmt-Moderately turbid (503) 
WAt-Highly turbid/very shallow (504) 

Wlm-Topographically low marsh (621) 
Whm-Topographically high marsh (622) 
Wtf-Tidal flat (623) 

Wga-Seagrasses and algal flats (624) 
Ws-Vegetated dredge spoil (625) 

B-Beaches (72) 

Bd-Dunes (731) 

Bds -Dredge spoil -barren (732) 

Bu -Undifferentiated barren land (77) 


Numbers refer to system in Anderson, et al . , 1976, U.S.G.S. Professional 
Paper 964. Two digit numbers identify Level II categories equivalent 
to U.S.G.S. Level II categories, although the descriptive terms may 
vary. Three digit numbers identify Level III categories. 
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3,4 Preliminary Evaluation of Techniques foir Analysis t>f LANDSAT 


Data 

Results of Interpreting LANDSAT Images of the Coastal Zone 
Indicate that a substantial degree of success can be achieved In 
mapping land cover/land use from LANDSAT film products, A knowledge 
of coastal geologic processes and biologic assemblages Is essential 
to this process because it enables the human interpreter to use 
much more than just reflectance in delineating coastal features. 

The shape of an object, its internal texture, and Its characteristic 
position with respect to adjoining environmental units can supercede 
reflectance as the basis for making classification decisions. Per~ 
haps because LANDSAT data originate In a digital format and are 
readily pnicessed by machine, these aspects of imagery interpre-' 
tatlon have received less attention than they certainly warrant. 

Reflectance alone, seen as the color tones of the false^color 
composite or the gray tones of a single-band Image, Is not absolute 
for identifying each type of land cover and land use. The growth 
phase of the vegetation, recent weather conditions, atmospheric 
conditions at the time of Image acquisition, and tide level are 
factors which must be taken Into account In Interpreting reflectance 
in a LANDSAT scene of the coastal region. The Interpreter using 
standard LANDSAT products therefore must rely on a familiarity with 
the coastal environments In order to compensate for the limit 
of resolution of about 80 m. and the 1:1,000,000 scale of the imagery. 
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Preliminary results of classifying LANDSAT tapes using un- 
supervised automatic data processing (ADP) techniques indicate 
that some classes of similar reflectance, such as beach, barren 
spoil and some urban/industrial, also can be discriminated only by a 
human interpreter using supportive Information. Advantages of 
ADP classification maps include 1) a large-scale display (1:24,000), 
when needed, 2) registration to U.S.G.S. 7 1/2 minute topographic 
maps with a RMS of 80-100 m., and 3) the use of statistical techniques 
to automatically reduce the volume of data to produce a map with 
15 to 20 spectral classes. 

3.5 Changes in ADP Techniques 

During this quarter, several ADP features have been added 
or modified to enhance the accuracy or speed of the classification 
process or display of the classified results. 

A very useful program, HGROUP, was acquired from the Regional 
Applications Project staff at NASA/JSC (Appendix B). When Classi- 
fying a LANDSAT scene for this investigation, Interest is usually In 
a small number of basic catagorles; e.g., water, wetlands, rangeland, 
urban, agricultural land, woodlands, and barren land. Even allowing 
subgroups of Interest within some of these basic categories, the 
result wou’^d rarely be more than about 15 classes on one map area. 
However, the use of ISOCLS clustering on a LANDSAT scene usually 
produces upwards of 25 to 35 subclasses. One thus needs a method 
of combining many of these subclasses to reduce the number of major 
classes. This often proves to be a difficult and tedious task. 


HGROUP is a valuable aid in this endeavor. Using a set of 
class means from ISOCLS that have been normalized, HGROUP performs 
a stepwise combining of classes. At each step, HGROUP combines 
the two classes that have the closest set of means. After com- 
bining two classes, a single set of means is generated for the com- 
bined class. The program continues combining classes until only 2 
classes remain. With the help of a cumulative error function 
generated at each step, HGROUP can be used to combine the original 
set of subclasses into the final set of classes of interest. 

Also, during this quarter, a program was written to extract 
boundaries from classified data (Appendix C). The boundaries are 
extracted in the form of chains of points which can be utilized by the 
Geographic Information System (GIS) now in development as part of 
the Texas Natural Resources Information System. Boundary extraction 
enables the G^S to produce pen plots at a smaller scale than 1:24,000 
from the LANDSAT boundary files, without the loss of information that 
occurs in line printer plots when lines and samples are dropped. At 
this time, polygons must be manually identified and labled, a time 
consuming task. Further capabilities (such as semi-automatic editing, 
polygon shading, area calculations, and overlaying of other base 
files) will become available as the GIS development progresses. 

Another program being written is one to help clean up "noise" 
in •*‘he classified data. This is being accomplished by changing 
the classification of a pixel to reflect that of its neighboring 
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pixels. The result should lead to more homogeneous areas in the 
final map. 

A program now in the design stage will use elements from 
the DAM package to perform simultaneous registration of two 
LANDSAT classification files of the same area (two different 
dates). As a matching pair of scan lines are being registered, 
they will be compared pixel by pixel for change. Only those 
pixels that have changed will be printed, thus giving a map 
displaying change between the two dates. 

After consulting with Regional Applications Project staff 
at NASA/JSC, the methods of picking control points has been 
altered. It was pointed out that picking all the control points 
within the topographic quadrangles to be classified would not 
produce registration for that area as accurately as using control 

points from the entire scene or at least from one entire CCT. 

The entire scene will be used in the future wherever possible. 

The LANDSAT ADP classification schedule was further re- 
fined to reflect addition of the HGROUP step and several minor 
changes in procedure v^hich have evolved during the work on test 
sites 2 and 5. This new schedule, shown in Table 2, will be 
used to support the cost analysis of work performed on site 4 by 
reference to the step number when recording time and equoiment 
utilization. 


17 



Table 2. 

LANDSAT Task 2 . Information Extraction: ADP Steps 

Scene ID: 

Test Site: 

Descri ption: 

1. Select LANDSAT scene and determine data tapes ID number. 

2. Examine available imagery. 

3. Merge data tapes or duplicate tapes If necessary. 

4. Estimate scan line and sample numbers for the areas of Interest. 

5. Generate grayscale maps of the area, (GRAYMAP/PICOUT) 

6. Obtain meteorological data 

7. Participate In orientation field trip. 

8. Establish control network (COEF) 

9. Classify water using DAM. 

10. Cluster all training areas within the scene (ISOCLS). 

11. Examine class statistics. 

12. Refine a training class If indicated by step 10. 

13. Use class statistics to build the look-up table (ELLTAB TABLE). 

14. Combine subclasses for display purposes (HGROUP). 

15. Classify the area (ELLTAB CLASSIFY). 

16. Register and display the classified results (REGISTER). 

17. Outline or color code homogeneous areas. 

18. Examine the classification map. 

19. Stop if satisfied with the results. 

20. Retrain on unclassified or poorly separated areas (ISOCLS). 

21. Go to step 12. 
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3.6 Summary of Accomplishments to Date 


As a consequence of the current LANDSAT project activities, 
certain results have been developed which appear significant. 

1) A land use/land cover classification scheme, that is 
representative of our investigation priorities, and 
applicable to the entire Texas Coast, has been developed 
which can be supported by LANDSAT interpretation. This 
scheme appears adequately detailed for use in coastal 
applications of several agencies. 

2) An update of the Bureau of Economic Geology's Environ- 
mental Geologic Atlas of the Texas Coastal Zone (Brown, 
project coordinator) has been completed, using (February 
1975 NASA photography (Mission 300), to provide a base map 

of the Coastal Zone features being studied with LANDSAT data. 

3) Systematic procedures have been identified and developed 
for the interpretation oF LANDSAT imagery according to 
the classification scheme. These image interpretation 
procedures can be used to generate a land use/land cover 
map of the coastal zone using conventional image inter- 
pretation and cartographic techniques, from which certain 
changes could be identified and evaluated for impact on 
the management of Coastal Public Lands. 



4‘) Automatic data processing (ADP) procedures have been im- 
plemented for the generation of computer classification 
maps, that are scaled and registered to 1 :24,000 U.S.G.S. 
topographic quadrangle maps. 

5) During this summer, information needs of the General 
Land Office in a specific coastal area will be evaluated 
with respect to the LANDSAT classification products 
developed. It is anticipated that LANDSAT products could 
play a significant supportive role in the GLO decision- 
making processes related to the management of the States' 
interest in the coastal public lands. 

In suRWjary, specific procedures have been developed for manual 
interpretation and computer processing of LANDSAT data according 
to a classification scheme which is relevant to Texas Coastal 
Public Lands management. Also, it appears feasible that an 
ongoing program, applying these procedures, could be established 
to produce a coastal land use/1 and cover map, on demand, to 
planning and management activities in the General Land Office and 
other Texas agencies. 

3,7 Program for Next Reporting Interval 

During the interval June-August 1976, a test of procedures de- 
veloped for image interpretation and ADP classification of LANDSAT 
data will be conducted in site 4, the Harbor Island area near Corpus 
Christi Bay, Texas, Four LANDSAT scenes will be analyzed for site 4 
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with the following dates: 2/25/75 (scene 1), Ifin^ (scene 2), 

7/10/75 (scene 3) and 12/16/72 (scene 4). Three scenes are winter 
scenes because cloud cover, obscuring part of test site 4, occurred on 
most scenes collected during the warmer seasons. The winter scenes 
will be used to Investigate changes in land use/land cover that might 
have resulted from construction along the bay margins and dredging 
of ship channels between 1972 and 1976. The summer scene will be 
used to compare differences in classification results between summer 
and winter scenes. 

A tentative schedule for this site 4 test Is shown in Table 3. 
Task codes listed for ADP and II refer to steps on the LANDSAT ADP 
Classification Schedule (Jones, et al., March 1976, Table 2, p. 22) 
and the Schedule for Image Interpretation Analysis ( Ibid , Table 3, 
p. 24). Estimates for completion of various tasks were provided by 
BEG and TWDB/TNRIS staff from experience gained from previous 
analyses for sites 2, 3, and 5, and are tentative depending on un- 
forseen delays. The last two tasks on the schedule will follow 
completion of the image interpretation and ADP analyses. In September 
and October. The analysis of problems/information of interest to 
the GLO will begin when at least two scenes are completed for ADP 
and image interpretation. Special interests for the GLO include 
the following topics: 

a) What is the area of inundation of tidal flats, marshes at 
different water stages? - such as water/1 and maps for each 
scene compared to one or more classification maps. 
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b) llhat is the "history” of ship channels as de~crifaed on 
LANDSAT scenes? - including turbidity plumes from dredges, 
growth of spoil islands, estimated acreage of bay bottom 
and marsh or grassflat covered, vegetation appearing on 
spoil, changing land use on adjacent shores. 

c) Have there been any changes in area or conditions of 
"wetlands" not related to dredging? seasonal changes? 

any trends" What impacts have the observed residential and 
industrial expansion had on Harbor Island area? Can 
expansion of these uses be related to the number of 
applications for easements, leases and structure re- 
gistrations in GLO? 

4.0 SIGNIFICANT RESULTS 

The image interpretation mapping techniques have been successfully 
applied to test site 5, an area with a semi-arid climate and consequently 
a very different coastal environment than the previously-studied test 
sites. As a result, the land cover/land use classification developed 
for this investigation required further modification. The system shown 
in Table 1 (Section 3.3.2) is applicable to the entire Texas coa?', and 
will be used during the remainder of this investigation. 

A new program, HGROUP, added to the ADP classification schedule 
(Table 2, Section 3.5) provides a convenient method for examing the spectral 
similarity between classes. This capability greatly simplifies the task 
of combining 25-30 unsupervised subclasses into about 15 major classes 
that approximately correspond to the land use/land cover classification 
scheme. 
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5.0 PUBLICATIONS 

A paper describing the techniques for mapping from standard 
LANDSAT film products, the land cover/1 and use classification and 
selected test site r>esults has been submitted to the Gulf Coast 
Section of the Society of Economic Paleontologists and Mineralogists. 
The paper will be published in October, 1976, at which time an oral 
presentation will be made. The content of the paper, which Is entitled 
"Interpretation of Unenhanced LANDSAT Imagery for Wetland and Land Use 
Delineation In the Texas Coastal Zone" by Rroert J. Finley, has been 
covered In this and previous quarterly progress reports. 

b.O RECOMMENDATIONS 
None 


7.0 FUNDS EXPENDED 


GENERAL LAND OFFICE (GLO) 

Labor 

Overhead 

Travel 

TOTAL EXPENDITURES REIMBURSED DURING THE 
4th QUARTER 

BUREAU OF ECONOMIC GEOLOGY (BEG) 

TOTAL EXPENDITURES REIMBURSED DURING THE 
4th QUARTER 

TEXAS PARKS AND WILDLIFE DEPARTMENT (TPWDl 

TOTAL EXPENDITURES REIMBURSED DURING THE 
4th QUARTER 


$4,606.00 

40.16 

160.00 

$4,806.16 


$ 


0 


I 0 
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TEXAS WATER DEVELOPMENT BOARD (TWDB) 

Labor $4,233.00 

Computer 600.00 

TOTAL EXPENDITURES REIMBURSED DURING THE 

4th QUARTER $4,833_.0_0_ 

CONSULTANTS 

Dr. Holz $2,250.00 

Dr. Schell 436.96 

TOTAL EXPENDITURES REIMBURSED DURING THE 

4th QUARTER $2,686.96 

CUMMULATIVE TOTAL EXPENDITURES REIM- 
BURSED DURING THE 4th QUARTER $12,326.12 

8.0 DATA USE AS OF MAY 31 , 1976 



IMAGERY 

Account #23790 
Amount 

CCT 

Account #G B3790 
Amount 

Value of Data 

$2,900.00 

$5,400.00 

Allowed 



Value Ordered 

> 

$1,482.00 

$1,800.00 

Value Received 

$1,460.00 

$1,800.00 


AIRCRAFT 
Account #G W3790 
Amount 


$9,588.00 

$9,564.00 

$9,564.00 



BALANCE 


$1 ,418.00 


$3,600.00 


$ 24.00 








9.0 AIRCRAFT DATA 

All NASA aircraft data flown for this investigation is now being 
used in a qualitative evaluation of map unit classification in sites 
2 , 3, and 5. The 1:40,000 color infrared photography of October 1975 
(NASA Mission 325) has been used extensively because of its good clarity 
and large scale. Some of the individual units evaluated with NASA 
photography in sites 2 and 5 are those described in the discussion of 
image interpretation results (Section 3.2). 
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APPENDIX A 


LANDSAT COVERAGE OF THE TEST SITES 2, 3, 4, 5 
FOR LANDSAT INVESTIGATION #23790 


ACQUISITION 

STATUS 

SCENE ID 

DATE 

CLOUD 

COVER 

QUALITY 

2 

Test Site 2: 

Sutniner: June - August 

1037 - 16251 

08/29/72 

20% 

8888 


1343 - 16253 

07/01/73 

20% 

8888 


1361 - 16252 

07/19/73 

20% 

8888 

2, 4 

1703 - 16175 

06/26/74 

10% 

8858 

2 

Fall: Sept. - Nov. 

1073 - 16251 

10/04/72 

30% 

8888 

2 

Winter: Dec. - Feb. 

1217 - 16261 

02/25/73 

20% 

8888 

2 

1505 - ■;S230 

12/10/73 

0P% 

2822 

2 

1901 - 16110 

01/10/75 

10% 

8808 

2, 4 

2375 - 16112 

02/01/76 

00% 


2 

1576 - 16152 

02/19/74 

00% 

8888 


Spring : Mar. - May 



1253 - 

16262 

04/02/73 

20% 

8888 

2, 4 

1289 - 

16261 

05/08/73 

00% 

8888 

2 

2051 - 

16140 

03/14/75 

00% 

8855 

2 

5027 - 

16050 

05/16/75 

10% 

5588 


ACQUISITION 

STATUS 


CLOUD 


\TUS 

SCENE ID 

DATE 

COVER 

QUALITY 


Test Site 3: 





Summer: 

June - Aug. 





1343 

- 16253 

07/01/73 

20% 

8888 


1361 

- 16252 

07/19/73 

20% 

8888 


1038 

- 16305 

08/30/72 

20% 

8888 


1362 

- 16305 

08/30/72 

20% 

8888 

2, 4 

1703 

- 16175 

06/26/74 

10% 

8858 


Fall: : 

Sept. - Nov. 





1092 

- 16312 

10/23/72 

20% 

8888 


1110 

- 16313 

11/10/72 

00% 

8888 


1452 

- 16291 

10/18/73 

00% 

7828 


Winter: 

Dec. - Feb. 




2, 4 

1146 

- 16314 

12/16/72 

00% 

8888 


1164 

- 16312 

01/03/73 

10% 

8888 


1182 

- 16313 

01/21/73 

00% 

8888 

2, 4 

2034 

- 16200 

02/25/75 

00% 

8888 


2016 

- 16200 

02/07/75 

10% 

5888 

2 

1578 

- 16264 

02/21/74 

10% 

8282 


Spri ng : 

Mar. - May 





1253 

- 16262 

04/02/73 

20% 

8888 


1289 

- 16261 

05/08/73 

00% 

8888 


1236 

- 16320 

03/16/73 

10% 

8888 


0 7- A-2 HiiPKODUClBlLlTY OF THE 

ORIQCNAL PAGE IS POOR 





ACQUISITION 

STATUS 

SCENE ID 

DATE, 

CLOUD 

COVER 

QUALITY 


1290 - 16315 

05/09/73 

00% 

8888 


1308-16314 

05/27/73 

20% 

8888 

2, 4 

1614 - 16261 

03/29/74 

10% 

8888 


1974 - 16133 

03/24/75 

00% 

8858 

2 

5028 - 16104 

05/17/75 

10% 

8885 

2 

Test Site 4: 

Summer: June - Aug. 

1326 - 16315 

06/14/73 

10% 

8888 

2 

1740 - 16225 

08/02/74 

20% 

8888 


1758 - 16221 

08/20/74 

20% 

8888 

2, 4 

5082 - 16080 

07/10/75 

10% 

8888 

2 

Fall: Sept. - Nov, 
1092 - 16314 

10/23/72 

10% 

8888 


1110 - 16320 

11/10/72 

10% 

8888 


1452 - 16293 

10/18/73 

10% 

8828 

2 

2268 - 16184 

10/17/75 

00% 

5555 

2, 4 

Winter; Dec. - Feb. 
1146 - 16320 

12/16/72 

20% 

8888 


1164 - 16315 

01/03/73 

20% 

8888 

2, 3* 

1182 - 16315 

01/21/73 

00% 

8888 


2016 - 16202 

02/07/75 

00% 

5885 

2, 4 

2034 - 16202 

02/25/75 

00% 

8888 

2, 4 

2376 - 16172 

02/02/76 
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ACQUISITION 

STATUS SCENE ID 

Spri ng : Mar. - May 

1236 - 16323 
1254 - 16323 
1290 - 16321 
5334 - 15523 
1308 - 16320 
1974 - 16135 
5028 - 16111 
Test Site 5: 

Summer : June - Aug. 

1362 - 16315 
1380 - 16314 
1722 - 16235 
1740 - 16231 
1758 - 16223 

Fal 1 : Sept. - Nov. 

1110 - 16322 
1776 - 16215 
1452 - 16300 

Winter : Dec. - Feb. 

1182 - 16322 
1506 - 16293 
2034 - 16205 


DATE 

CLOUD 

COVER 

QUALITY 

03/16/73 

20% 

8888 

04/03/73 

10% 

8888 

05/09/73 

20% 

8888 

03/18/76 



05/27/73 

10% 

8888 

03/24/75 

10% 

8858 

05/17/75 

10% 

5588 


07/20/73 

20% 

8888 

08/07/73 

20% 

8888 

07/15/74 

20% 

8888 

08/02/74 

10% 

8888 

08/20/74 

10% 

8888 


11/10/72 

10% 

8888 

09/07/74 

20% 

5855 

10/18/73 

20% 

8888 


01/21/73 

00% 

8883 

12/11/73 

10% 

8888 

02/25/75 

00% 

8888 


1 

2 

2 


2, 4 
2 

2 

2 


2. 4 


2, 4 
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ACQUISITION 

STATUS 


SCENE ID 


DATE 


CLOUD 

COVER QUALITY 


Spring : Mar. - May 


1614 - 

16270 

03/29/74 

20% 

8888 

1974 - 

16142 

03/24/75 

10% 

8888 

2070 - 

16203 

04/02/75 

20% 

8588 

1290 - 

16324 

05/09/73 

20% 

8888 


ACQUISITION STATUS 

1 Imagery on Order 

2 Imagery on Hand 

3 Tapes on Order 

4 Tapes on Hand 

* Tape could not be reproduced by EDC/Goddard 
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DB02Q0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 
3d 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


-02*ELLTAB.HGROUP 


C 

C 

c 

c 

c 

C‘ 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM HGROUP 

HIERARCHICAL PROFILE-GROUPING ANALYSIS* 

PARAMETER CONTROL-CARD FIELDS. 

COL 1-5. number of VARIABLES (MAX = 54). 

COL 6-10. NUMRER OF SUBJECTS {MAX = 54). 

COL U-15. LEVEL OF GROUPING TO BEGIN GROUP-MEMBERSHIP PRINTING. 
COL 20. 1 = STANDARDIZE DATA ON EACH VARIABLE BEFORE GROUPING. 
COL 25. 1 = TRANSPOSE DATA MATRIX IN ORDER TO GROUP VARIABLES. 
FORMAT MUST SPECIFY AN ALPHANUMERIC SUBJECT-CODE FIELD » FOLLOWED BY 
NV SCORE FIELDS. IF DATA MATRIX IS TRANSPOSED (COL 25 = 1) # 
group-membership CODES WILL BE SERIAL NUMBERS OF VARIABLES. 
SUBPROGRAMS REQUIRED ARE SUMF AND CCDS. 

DIMENSION D{54rS4)# KG{54) ^ W(54)f KF(20) 

REAL+8 LC(54)f KCC54) 

ND = 54 

5 CALL CCDS (KFr NVf NSf KPf KSr KT) 

T = NS 

C READ ALL DATA CARDS AND STANDARDIZE COLUMNS (VARIABLES)# IF 
C OPTIONED 

DO 10 I::lfNS 

ID READ KFr KCil ) p (D(IfJ)r JslfNV) 

IF (KS .eg. 0) GO TO 20 
DO 15 JrliNV 

A = SUMFtDf Jr NSr ND) / T 
S = SQRTtSUMFtO# Jr -NSr NO) / T - A *1^ A) 

DO 15 I=lrNS 

15 D(IrJ) = (D(lrJ) - A) / S 
20 IF (KT .eg. 0) GO TO 30 
C TRANSPOSE DA'' A MATRIX# IF OPTIONED. 

N = MAXOCNSr NV) 

DO 25 I=lrN 
DO 25 J=lrN 
X = DCIrJ) 

O(IfJ) = D(Jrl) 

23 D(Jrl) = X 
NS = NV 
NV = T 

C CONVERT DATA MATRIX TO INITIAL MATRIX OF ERROR POTENTIALS. 

30 DO 45 I::lrNS 
DO 35 J=lrNV 
35 W(J) := D(IrJ) 

DO 45 J-IfNS 
DilpJ) - 0.0 
DO 40 KrlfNV 
40 D(I#J) = D£IrJ) 

45 D(IrJ) = D(IrJ) 

DO 55 IslrNS 
DO 55 J=IfNS 
55 D(Jrl) r 0,0 


(D(JrK) - WCK))**2 

2.0 


52 NG=N5 

53 C INITIALIZE GROUp-MEMBERSHIP AND GROUP-N VECTORS. 

54 DO 60 I::lrNS 

55 KGfDrl 

56 63 Wm=1.0 


^ 7 - B-1 


57 

58 

59 

60 
61 
6a 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 
75 

74 

75 

76 

77 

78 

79 

80 
81 
82 
83 
64 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 
113 


C LOCATE optimal COMBINATIONf IF MORE THAN 2 GROUPS REMAIN. 

65 NG=NG-1 

IF (NG #EQ* 1) GO TO S 

x=io.o*no 

DO 75 121 f NS 

IF (KG(I> .NE. I) GO TO 75 
DO 70 J=IfNS 

IF Cl .EQ. J .OR. KG(J) .ME. J) GO TO 70 
DX 2 DCIfJ) - DCIfI) - D(J#J) 

IF (DX .GE. X) GO TO 70 

X2DX 

L=I 

M=J 

70 CONTINUE 
75 CONTINUE 
NL = W(U 
NM = W(M) 

WRITE (6p 80) NGf Lt NLf NMf X 
flOOFORMAT (/ 14 r 25H GROUPS AFTER COMBINING Gf I5f 

14H CN=f I3f 7H) AND Gr I3 f 4H (Nr# I3r lOH ) f ERROR -f 
2 F16.6) 

C MODIFY GROUP-MEMBERSHIP AND GROUP-N VECTORSi AND ERROR 
C POFENIALS. 

WS = W(L) + WCM) 

X 2 DCLfMI * WS 

Y 2 OCLfL) + WCL) + DCMfM) ♦ W(M) 

OCLiU 2 O(LfM) 

DO 85 I21,NS 

IF (KGCn .eg. M) KG(I) 2 L 
85 CONTINUE 

00 95 121 F NS 

IF (I .EQ. L .OR. KG(I) .NE. I) GO TO 95 
IF £I .6T. L) GO TO 90 

ODCIfL) = £D(IiU * CW(I) + W(D) + D(IfM) tWCl) + WCM)) 

1+ X - Y - OCIfI) * WCD) / (WCI) + WS) 

GO TO 95 

90DD(LfI) 2 CDCLfI) * CW(L) ^ W(I)) ^ (DCHfI) + D(IfM)) 

(W(M) + W(I)I + X - Y - DCIfI) W(D) / tWCl) + WS) 

95 CONTINUE 
W ( L) ~ WS 

IF (NG .GT, KP) GO TO 65 

C PRINT group memberships OF ALL OBJECTS* IF OPTIONED. 

DO 115 I=1*NS 

IF (KG(I) ,NE. I> GO TO 115 
L=G 

00 100 J=I*NS 

IF (KG(J3 .NE. I) GO TO 100 
L = L + 1 
LC(L) = KCtJ) 

IF (KT .EG. 1) LC(U = J 
100 CONTINUE 

IF CKT .EQ. 0) GO TO 102 
IF (KT .EQ. 1) GO TO 104 

104 WRITE (6il05) I* L* (LC(J)* J=1»U 

105 FORMAT (2H G, 13* 4H (N=* 13* 2H) » 2514 / (14X* 2514)) 
102 WRITE (6*110) I* L* (LC(J)* U=1*U 

no format (2H G* 13* 4H (N=* 13* 2H) » 15A7 / 114X* 15A7)) 


114 

115 

116 


115 CONTINUE 
GO TO 65 
END 
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DB0200-02*ELLTAB. NORMAL 

1 c***program normal 

2 C THIS PROGRAM NORMALIZES A VECTOR ARRAY 

5 C 

4 C NV = VECTOR DIMENSION 

5 C NS = NUMBER OF SAMPLES 

6 C ■ KF s: VARIABLE FORMAT FOR READING - CLASS NAME AND DATA VECTOR 

7 C 

8 DIMENSION KF<2fl ) ,KC (50) ^D(50 1 50) f VMAXC50) 

9 C READ PARAMETERS # FORMAT AND DATA ARRAY 


10 



READ(5flOO) NVtNSrKF 

11 


100 

F0RMAT(2I5*/i-20Ait) 

12 



DO 10 T=1»NS 

13 


10 

READ(5fKF) KC C I) , tD<I » J) f J"! f NV) 

14 

C 

find maximum value for each variable 

15 



DO 20 J=lfHV 

16 



DO 20 I=1»NS 

17 


20 

VMAX(J)=AMAXl(VMAX(U)tD(ItJ) ) 

18 

C 

NORMALIZE THE VECTORS 

19 



DO 30 J=1»NV 

20 



DO 30 1=1 tNS 

21 


30 

0(1»J)=D(I»J)/VMAX(J) 

22 



DO 40 l=ltNS 

23 


40 

WRITE (1*KF) KCm. (D(I»J) »Jsl»NV) 

24 



WR1TE<6»200) 

25 


200 

FORMAT(//i* normalization COMPLETED’ »//) 

26 



STOP 

27 



END 

QBRKPT 

PRINTS 




DB0200-02*ELLTAQ . SUMF 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTlOi SUMF (X^ KK r MNr ND) 

COMPUTES SUM X OR SUM X:*:+2 FROM A VECTOR. 

X ARRAY CONTAINING THE SCORES TO BE USED, 

KK = ROW OR COLUMN NUMBER IF X IS A MATRIX. SET = 1 IF X IS A VECTOR. 
IF KK IS positive AND NOT IT IS A COLUMN VECTOR. 

IF KK IS negative AND NOT 1 » IT IS A ROW VECTOR. 

NN = NUMBER OF VALUES TO BE SUMMED. IF NEGATIVE^ SUM X+^2 COMPUTED. 

ND r: NUMBER OF ROWS (OR ELEMENTS) DIMENSIONED FOR X IN THE 
CALLING PROGRAM. 

DIMENSION X(54rl) 

SUMF =0.0 
N = lABStNN) 

K = lABS(KK) 

IF (NN? 5*55^10 
5 IF (KK) 15fS5f25 
10 IF (KK) 35r55»45 
15 DO 20 I=lrN 
20 SUMF = SUMF + X£K^I)*^2 


21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 


RETURN 

25 DO 30 I=lrN 
30 SUMF = SUMF + X(I?K)*>»c2 
RETURN 

35 DO 40 I=lrN 
40 SUMF = SUMF + XCKfl) 
RETURN 

45 DO 50 I=lrN 

50 SUMF = SUMF + X(1#K) 

55 RETURN 
END 
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APPENDIX C 
PROGRAM EXTRACT 




V.v 
\ ■■ 
K 




h 



■S' 


f 

I 

\ 

i 


;?7-c. 


J^^0200*UAM[)SAT. EXTRACT 


1 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
30 
59 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 
5to 


C***PROGRAM EXTt>ACT 
C 

C THIS PROGRA:^ extracts LiOUUUARlES FROM LANDSAT CLASSIFICATION FILES 
C UNIT 8 = registered MAP USED AS INPUT 
C UNIT 15 = OUTPUT IGF BOUNDARY MAP 


dimension LIN^^l (4000 )rLINE2 14000) i LPRTl (4000 ) fLPpT2 (4000 ) 
DATA IN#IOlJT/5»15/aDrNPTS#XINC/lflf0.1/tISAVF/0/ 
YINCr-l*/6* 

CALL SETAOPdOUT,!) 

C CALL SETMRG(2MM»66rO0r00 . O 
C**+READ IN CONTROL POINTS 

5 READ(5r5OOrENO=6fERR=05) NROW»NCOL 
500 FORMAT () 

Xl-tNCOL-l)i^XlNC + XlMC/2. 

Y3 = Cr!ROW“l)*YlMC + YlNC/2. 

WRITE(IOUT) IDr.-JpTSfXl»Yl 

ID-ID+I 

GO TO 5 

6 NPTS=2 

C**^INITIAL12E ARRAYS 
00 10 I=if40gn 
LPRT2 ( I) =6H 
10 LINE2(I)=6H 

C***READ FIRST CLASSIFIED LINE INTO CORE 

READ(INfEN0=99fERR=86) hPTX » ( L lNE2 ( I ) rNPlX ) 

NROW=l 

IENO-MPIX-1 

C=M«tLOCATE ROUNDARY POINTS TN THE FIRST LlNF 
DO 20 T=1»IENT 

IF {LINE2(I)-LTnE2(I + l}) 14t20d4 
C***EXTRACT BOUNDARIES 
14 X1=I*XXNC 

Y1=(NR0W-1)+YTNC 

X2=I*XINC 

Y2-NROW*YINC 

LPRT2(I)riLlNL2n) 

LPPT2(I + l)rLIN^2n + l) 

’.vRlTF (TOUT) npNPTSrXlrYl^Xi»Y2 
I0=ID+1 
20 CONTINUE 

C***5HIFT DATAr THEN PROCESS ‘^EYT LINE 
30 DO 40 I=ir400n 
LPWTim=LPRT?f I) 

LPPT2(IJ=6H 
LTNEltI)=LlNE?( n 
40 LINE2(I)=6H 

REAO( rJrENn=9a»ERR=0ol iPIXi (L INE2 ( r 1 = 1 ^ MPI X) 
r*DOW=NRO^/^ + l 

C=^*tL0CATE ^»OUNDARY ^r.'lNTS 
DO 50 

IF (IpEl^*NPIX) <yQ TO tift 
IF (LTNE2(T)-LTNL2U + 1) ) 44»u#i.r4*+ 

C***EXToaCT r^o'J^inAPTr-^ 

44 Pl = I*XTr,C 
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57 

58 

59 

60 
61 
62 

63 

64 

65 

66 
67 
6B 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
83 
04 

85 

86 
87 
85 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 
qg 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 
113 


P2=I*)(INC 

Q2rNROW*YlNC 

LPRT2nJ=rLINb?(I) 

LPRT2 ( I+l ) =LIHF2 (1+1) 

WRITEdOUT) lOrNPTSrPl PtiltP2r02 
ID=ID+1 

C*^*CHECK FOR BOUNDARY POINTS WITH PREVIOUS LINE 
46 IF (UrME2(I)-LTNEl{I)) 48r50f48 
48 LPRTKIJ-LINEKI) 

LPRT2(IJ=:LINE2(I) 

IF (ISAVE^EQ.O) GO TO 150 

C*+«SAVE NEW boundary LINE 
XXl=(I-l)^>flNC 
YY1=(NR0W-1)*YTNC 
XX2=1*XINC 
YY25!(NR0W-1)*YTUC 

C*i>*IF THE END POINTS DON'T MATCH WRITE THE CHAIN 
IF (X2,E0.XXl.AnD.Y2*FQ*YY2) GO TO 145 
WRITE (TOUT ) in # r XI f Y1 » X2 r Y2 
ID=ID+1 
X1=XX1 
Yi=YYl 
X2=!XX2 
Y2=YY2 
60 TO 50 

C***TIE ADJACENT BOU^'OAPY CMAINS TOGETHER 

145 X2=XX2 
Y2::YY2 
GO TO 50 

A NErt CHaT^' 

150 X1 = {I-1)1^XINC 

Yl=(WOW-l)*YrMC 
X2rI#XINC 
Y2=(NROW“l)*YtNC 
ISAVE.^1 
50 CONTINUE 

C»»*OUTPUT BOUNDARY ON LINE PRlNTERr THEN GET HEVT LiNr 
CALL MAP 

IF (ISAVE.EO.ni on TC 3'i 

^^PITE(IOUT) V>f -PTS»XUY1#X2>Y2 

IO=ID+l 

ISAVE=0 

GO TO 30 

C***KRROR READING CO^^KOt POINTS 

85 WRITE(6»101) 

ICJl FORMATC* ERROR HEAPING CONTROL POINTS') 

STOP 

C***ERROH READING FILE 

86 wRITr(6il00) 

ion FORMATC EPRO'^ PEAPirjrt CLASSIFICATIOi. FILE') 

STOP 

C*»*END OF FILE 
99 ENHFILE IOMT 
CALL '-’AP 
WRlTEC6a99) TP 

199 FORMAT(///i IlOf » CHATS EXTRACTED*) 
v.RITE(6f?nn) 
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114 

U5 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 


200 FORMAT(/f» BO'JMUARY EIXTRACTION COMPUTED^ 
C call SETMF?G(2f »-*i66iD6r(J3 . M 
STOP 

C OUTPUT boundary *^AP ON LINE PRINTER 
SUBROUTINE 
NLINE=NLINE+1 
LA5T=MINO(NPIV»120) 

vJRlTE (31^300) MLlNEr (LPRTl { I ) ^ 1=1 1 LAST) 

300 F0RMAT(1X» J4rlYrl20Al) 

IF tNPlX.LT*12t) GO TO 10 
LAST=MIN0(MPIX»240) 

WRITE (32# 300} MLlNEi (LPRTl ( I) # 1=121 # LAST) 
IF CNPIX*LT,2413 GO TO 10 
LA5T=MlNO(NPlXr360) 

WRITE (33# 300) MLlNE# (LPKTK I) #1=241# LAST) 
IF (rJPIX,LT.361) GO TO 10 
LAST=VIN0 (NPiyr480) 

WRITE (34# 300) NLlME# (LPuTl (I) r 1=361 rLAST) 
IF (NPIX*LT.48l) GO TO 10 
LAST=*tIND(NPIX#600) 

WRITE (35# 300) NLINE# (LPRTl C I )# 1=481 # LAST) 
IF (NPIX*LT,601) GO TO 10 
LAST=‘^IN0(NPIXr720) 

10 RETURN 
END 
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